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Abstract 

A one-equation linear turbulence model and 
a two-equation nonlinear explicit algebraic stress 
model (EASM) are applied to the flow over a mul- 
tielement airfoil. The effect of the K-e and K-u> 
forms of the two-equation model are explored, and 
the K-e form is shown to be deficient in the wall- 
bounded regions of adverse pressure gradient flows. 
A new K-u> form of EASM is introduced. Nonlinear 
terms present in EASM are shown to improve pre- 
dictions of turbulent shear stress behind the trailing 
edge of the main element and near midflap. Cur- 
vature corrections are applied to both the one- and 
two-equation turbulence models and yield only rela- 
tively small local differences in the flap region, where 
the flow field undergoes the greatest curvature. Pre- 
dictions of maximum lift are essentially unaffected 
by the turbulence model variations studied. 

1 Introduction 

In recent studies, 1,2 analyses of high-lift multi- 
element airfoil configurations were performed to as- 
sess the predictive capability of turbulence models. 
These studies were undertaken in part to determine 
why, in comparison with nominally two-dimensional 
(2-D) experiments, CFD overpredicts maximum lift 
and the angle at which maximum lift occurs. In 
general, with the full airfoil configuration, it is diffi- 
cult to isolate the deficiencies of turbulence models. 
However, the studies revealed three areas where tur- 
bulence model predictions were possibly deficient and 
could adversely affect the overall prediction of the 
flow field. These areas w'ere (1) prediction of transi- 
tion location, (2) prediction of downstream evolution 
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of the slat wake, and (3) accounting for streamwise 
curvature effects in the flap region. Additional stud- 
ies over the last several years have examined some of 
these possible modeling deficiencies in isolation by 
using “unit problems.” Section 2 highlights results 
from these studies. 

The current study was an effort to gauge the ef- 
fects of recent turbulence model advances (some of 
which resulted from unit problem investigations) on 
a multielement airfoil near maximum lift. The study 
investigated three major areas: the effects of K-e 
vs. K-u formulations, the effects of nonlinear terms, 
and the effects of curvature corrections. Following 
section 3, in which the numerical method and turbu- 
lence models are briefly described, the three major 
areas are discussed in turn in section 4. For all re- 
sults in this paper, we have focused our attention on 
the region of the configuration near the main element 
trailing edge and flap, where previous studies have 
shown the greatest discrepancies with experimental 
data. 

2 Summary of Unit Problem Investigations 

Rumsey et al. 3 investigated the ability of three 
turbulence models to model the effects of convex 
curvature in a U-duct unit problem. In particular, 
convex curvature is known to cause a suppression of 
the turbulent shear stress in the outer part of the 
boundary layer. The three turbulence models — 
one-equation Spalart-Allmaras (SA), 4 two-equation 
Menter shear-stress transport (SST), 5 and two- 
equation explicit algebraic stress model (EASM) 6 
— all behaved similarly in the curved region, and 
all failed to predict the suppression of the turbulent 
shear stress. It was shown that a Reynolds stress 
model (RSM) can predict the suppression, and that 
the source of the error in EASM (which is derived 
directly from RSM) is the assumption of anisotropy 
equilibrium in the Cartesian frame of reference: 
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of the slat wake, and (3) accounting for streamwise 
curvature effects in the flap region. Additional stud- 
ies over the last several years have examined some of 
these possible modeling deficiencies in isolation by 
using “unit problems.” Section 2 highlights results 
from these studies. 

The current study was an effort to gauge the ef- 
fects of recent turbulence model advances (some of 
which resulted from unit problem investigations) on 
a multielement airfoil near maximum lift. The study 
investigated three major areas: the effects of K-e 
vs. K-lo formulations, the effects of nonlinear terms, 
and the effects of curvature corrections. Following 
section 3, in which the numerical method and turbu- 
lence models are briefly described, the three major 
areas are discussed in turn in section 4. For all re- 
sults in this paper, we have focused our attention on 
the region of the configuration near the main element 
trailing edge and flap, where previous studies have 
shown the greatest discrepancies with experimental 
data. 


2 Summary of Unit Problem Investigations 


Rumsey et al. 3 investigated the ability of three 
turbulence models to model the effects of convex 
curvature in a U-duct unit problem. In particular, 
convex curvature is known to cause a suppression of 
the turbulent shear stress in the outer part of the 
boundary layer. The three turbulence models — 
one-equation Spalart-Allmaras (SA), 4 two-equation 
Menter shear-stress transport (SST), 5 and two- 
equation explicit algebraic stress model (EASM) 6 
— all behaved similarly in the curved region, and 
all failed to predict the suppression of the turbulent 
shear stress. It was shown that a Reynolds stress 
model (RSM) can predict the suppression, and that 
the source of the error in EASM (which is derived 
directly from RSM) is the assumption of anisotropy 
equilibrium in the Cartesian frame of reference: 
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where 6,j = [r,j/(2/\ )] — (%/ 3) and where A' = 
r nn /2 is the turbulent kinetic energy. 

Gatski and Rumsey 7 showed that by assuming 
Eq. (1) to hold, not in the Cartesian frame but rather 
in the frame defined by the principal axes of the 
strain rate tensor, a new form of the EASM could 
be derived that takes into account the flow field cur- 
vature. With the new EASM curvature-corrected 
(EASMCC) method, the suppression of the turbu- 
lent shear stress near the convex bend of the U-duct 
was accurately predicted. 8 

A curvature correction for the SA model has also 
been developed by Spalart and Shur 9 and applied 
to the same U-duct flow in Shur et al. 10 This cor- 
rection, Spalart- Allmaras for rotation and curvature 
(SARC) is based similarly on the rate of change of 
the principal axes of the strain rate tensor, but it 
also includes a heuristic function f r \ (that multi- 
plies the model’s production term) not present in 
the EASMCC. 

A different unit problem of wake development in 
various pressure gradients (an experiment conducted 
by Liu et al. 11 ) was investigated by Carlson et al. 12 
This unit problem was motivated by the fact that 
most turbulence models have overpredicted the slat 
wake depth and width on the multielement airfoil, 
even when transition is specified according to ex- 
perimental measurements. 13 Using a linear two- 
equation eddy viscosity model and EASM, Carlson 
et al. showed that the isolated effects of pressure 
gradient on near-wake mean flow development could 
be reasonably predicted by both turbulence models; 
results compared favorably with the experiment. 

However, the pressure gradients experienced by 
the slat wake in the real multielement configura- 
tion are highly variable (a short, very strong favor- 
able gradient followed by a longer adverse gradient), 
and these gradients can be significantly stronger 
than those imposed in the unit problem experiment 
(which had constant dp/dx). Therefore, it is still 
not known whether (1) the stronger (or variable) 
pressure gradients in the multielement wake cause 
a failure of the turbulence models, or (2) the turbu- 
lence models are not to blame and some other effects 
not being modeled in the CFD, such as unsteadiness 
or three-dimensional (3-D) effects, are causing the 
discrepancy. 

Finally, on-going separate work is focusing on 
transition prediction within the context of the 
EASM, based on rigorous mathematical treat- 
ment. 14,15 However, this effort is quite complex 
and longer term. A shorter term, empirically 
based method for predicting transition within a two- 


equation turbulence model has been applied to the 
multielement configuration by Czerwiec et al. 16 

3 Numerical Method and Turbulence Models 

The CFD code used in the current investigation 
was CFL3D, 17 a widely used structured-grid upwind 
finite- volume method. Details about the code can be 
found in the user’s manual referenced. 

Of the three fundamental investigations described 
in this paper, the first is the effects of the K-e 
vs. K-u) underlying formulation for the EASM tur- 
bulence model. Two versions of the EASM are 
denoted by EASM(A'-e) and EASM(A-cj), respec- 
tively. The EASM(AT-cj) is described in detail in 
the appendix, whereas the EASM(A'-£) has been de- 
scribed previously in Ref. 3 and will not be repeated 
here. However, note that the EASM(A'-e) has an 
additional minor modification, described by Gatski 
and Rumsey, 7 in order to correct for a tendency of 
the original model to produce excessive levels of eddy 
viscosity near the center of wakes; these excessive 
levels result in a nonphysical, local “flattening” of 
the velocity profiles. Equation (4) in Ref. 3 has been 
replaced by 

9 = 7oj + 7i* > (2) 

where 

7o = 7o — 1 (3) 

and 

1,1_T, + 1 + l c.,-1 )■ 

This modification only slightly affects results for the 
log layer, where C ^ (Eq. (29) of Ref. 3) now re- 
quires the value of 0.0885 for EASM(/i-£). The 
EASM(AT-w) described in the appendix also includes 
the above modifications to g (= cti/r). 

The second fundamental investigation described 
is the effects of the nonlinear terms themselves. In 
EASM, the turbulent stress tensor is not only a func- 
tion of the strain rate, but of two additional nonlin- 
ear terms as well (see the appendix) . For this part of 
the study, we used the EASM(A'-u;) with and with- 
out its nonlinear terms in place. 

The third fundamental investigation is an explo- 
ration of the effects of curvature on the multiele- 
ment airfoil case considered. Four turbulence mod- 
els were employed for this part of the study: SA, 
SARC, EASM(A'-u;), and EASMCC(A'-w) (in other 
words, SA and EASM(A-w) with and without their 
respective curvature corrections). 

Brief descriptions of the curvature corrections for 
SARC and EASMCC are given here. Both use the 
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kinematic strain rate and rotation rate tensors, Sij 
and W{j respectively, defined in the appendix. These 
corrections are similar in that they both employ the 
Lagrangian derivative of the strain-tensor principal 
axes, given by 


Da _ 1 

Dt ~ 2(S? 1 + S? a ) 



DSu 

Dt 


~ S \2 


DSn " 

Dt 


( 5 ) 


However, the two curvature corrections were devel- 
oped independently and are different in many other 
respects. 

In the current implementation of the SARC 
model, a portion of the SA model’s production term, 
Cm [1 - ft<i\Wv, was replaced by cm[/m - fti] Wv, 
where 


2 r* 

fn = (1 + c r i) [l - c r 3 tan _ 1 (c r 2 r)] - c rX , 

. W 

and Cri = 1, c r 2 = 12. The constant c r 3 was assigned 
to be both 1.0 and 0.6 in Spalart and Shur , 9 who ad- 
mitted they are still experimenting with the heuris- 
tic function f r \. For the current study, we used 
c r3 = 0.6. The function r* is given by r* — S/W } 
where S = ^SijSij and W = y/ 2 W~JW~. For 2-D 
flows and no system rotation, the expression for r in 
Ref. 10 can be reduced to 


(Da\ &W 12 (Sh + S 2 12 ) 
\ Dt ) D A 


(?) 


where Da/Dt is given by Eq. (5) and D = 

v/0.5 (S 2 +V 2 ). 

In the study of non-Newtonian constitutive rela- 
tions (e.g., Schunk and Scriven , 18 Souza Mendes et 
al. 19 ), a measure of relative rotation rate is based 
on the principal axes of the strain rate tensor. As 
mentioned in the introduction, Gatski and Rumsey 7 
used this measure to derive EASMCC, which takes 
into account the flow field curvature. 

The method for implementation of EASMCC in 
2-D is as follows. The rotation rate tensor W t j in 
the model is replaced by Wij — where the 

constant 02 is defined by the pressure-strain corre- 
lation model (recall that Wij — Qij is the absolute 
rotation rate tensor). For the SSG model 20 used 
here, — (2 — C 4 )/ 2 and C\ — 0.4. The tensor Q*j 
is given by 


flij = 


0 Da/Dt 

— Da/Dt 0 


( 8 ) 


and Da/Dt is given by Eq. (5). 


4 Results 

In the current study we focused primarily on the 
region of the flow field in the vicinity from the main 
element trailing edge to the midflap area. This re- 
gion was identified by Ying et al . 2 as an area where 
CFD results generally differ substantially from ex- 
perimental data. It is also the region where any 
streamwise curvature effects would be expected, be- 
cause the flow turns rapidly through 30-40 deg 
as it passes over the flap. A diagram showing 
the 30P-30N multielement configuration is shown in 
Fig. 1, with the current region of focus delineated. 

Unless otherwise noted, all computations shown 
were performed on a “free- air” grid at a = 19 deg, 
M — 0.2, and Re = 9 million. The effect of modeling 
the lower and upper tunnel walls in the CFD grid 
was explored in Ref. 21 and is not repeated here. 
Also, grid effects have been studied previously for 
this configuration (see discussion on numerical sen- 
sitivity in Rumsey et al. 1 ); the current 4-zone grid 
contains 135,428 grid points. The transition loca- 
tions on each element were specified at the end of 
the measured ranges given in Ref. 13. 

4.1 Effects of 6 vs. a? Formulations 

The EASM model can be coupled with any two- 
equation model formulation; in this study, we ex- 
plored the effects of K-e vs. K-uj. Velocity profiles 
are shown in Fig. 2 at x/c — 0.85, near the trail- 
ing edge of the main element upper surface. In the 
plot, d refers to the normal distance from the airfoil 
surface and c is the stowed chord length. Results 
from using EASM(A"-u>) showed better agreement 
with the main element boundary layer thickness and 
slat wake depth than results with EASM (K-e). Al- 
though not shown, the latter model overpredicted 
the turbulence levels in the wall-bounded adverse 
pressure gradient regions over both the slat and the 
main element. This overprediction caused thicker 
predicted boundary layers, and consequently wakes 
too wide and deep. 

The inability of K-e models in general to handle 
wall-bounded adverse pressure gradients is a known 
problem characteristic of the 5 -equation. See, for ex- 
ample, Wilcox , 22 Rodi and Scheuerer , 23 and Nagano 
and Tagawa . 24 In Ref. 23, the shortcomings are ex- 
amined and it is shown that the generation term of 
the 5 -equation has to be increased. 

Figure 3 shows a plot of the predicted veloc- 
ities in wall variables as compared to Spalding 
theory 25 at x/c = 0.85. The EASM(A'-5) did 
not obtain the correct slope of the log layer in 
this adverse pressure gradient flow. (Although not 
shown, EASM(AT-5), like other K-e models, per- 
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forms well for wall-bounded, zero-pressure-gradient 
flows.) EASM (/\-w) had no trouble handling the 
adverse-pressure-gradient flow, and obtained the 
correct log-layer slope and position. 

From these results, it appears clear that 
EASM(A'-e) should not be used for flow fields of 
the type explored in this study. (For that mat- 
ter, any K-e model without a modification, e.g. 
Ref. 23, probably cannot correctly predict wall- 
bounded adverse pressure gradients flows.) There- 
fore, for all results in the remainder of the paper, 
only the EASM(A'-u>) version of EASM was used. 
For brevity, EASM(R'-w) and EASMCC(A'-uj) mod- 
els will henceforth be referred to as EASM and 
EASMCC, respectively. 

4.2 Effects of Nonlinear Terms 

Figure 4 shows experimentally measured 26 turbu- 
lent shear stresses at xjc = 0.85 near the trailing 
edge of the main element and at xjc — 0.898 on the 
flap. In the figure, the latter curve has been shifted 
to align the relevant features. The experiment indi- 
cated almost no change in the minimum u'v' level 
(near d/c = 0.01) between these two nearby loca- 
tions. On the other hand, the SA model predicted a 
large change in the turbulent shear stress as the flow 
passed from the trailing edge to the wake, as shown 
in Fig. 5. Ying et al. 2 surmised that the reason 
for this too-rapid change was due to the Boussinesq 
eddy- viscosity assumption inherent in the SA (and 
other linear models). In the Boussinesq assumption, 
the turbulent shear stress is assumed to be directly 
proportional to the strain, so that = —2 u t Sij. 
In this flow field, as the strain rate changes rapidly 
from the trailing edge to the wake, so too does the 
turbulent shear stress because v t from the SA varies 
smoothly and gradually there. 

In contrast, the EASM (Fig. 6) did a better job 
maintaining the minimum shear stress level between 
these two stations. But what was the role of the 
nonlinear terms in the improved predictions? This 
question was investigated by solving EASM as a lin- 
ear eddy viscosity model (i.e., by forcing the non- 
linear terms to be zero). Results, shown in Fig. 7, 
indicate a large difference from the results of Fig. 6. 
Thus, the nonlinear terms appear to be necessary to 
capture more closely the streamwise development of 
the turbulent shear stress in this region. 

4.3 Effects of Curvature Correction 

The curvature correction for EASMCC encoun- 
tered some numerical difficulties for this case, asso- 
ciated primarily with the explicit one-to-one inter- 
faces in the wake-cut regions of the C-grids. Near the 
connectivity interfaces, small discontinuities could 


result in large levels of DSij/Dt , which would feed 
back into the solution through and worsen the 
discontinuities. To remove this problem, the cur- 
vature terms were turned on only in regions of the 
grid that contain a wall at the k = 1 index loca- 
tion. Furthermore, the velocities used in the deter- 
mination of DSij/Dt were smoothed by using 10 it- 
erations of an explicit point Gauss-Seidel Laplacian 
smoother, 27 and the elements of the Qij matrix were 
limited to prevent unreasonably large levels from oc- 
curring. Note that these numerical difficulties did 
not occur with EASMCC for simpler (e.g., single- 
zone grid) cases. 7 Because of the smooth limiting 
inherent in the tan" 1 function in the empirical f r \ 
expression, SARC did not encounter any numerical 
difficulties for this case. 

Velocity profiles from the four turbulence models 
are shown at three stations, from the trailing edge 
of the main element to the middle of the flap, in 
Figs. 8-10. The corresponding turbulent shear stress 
profiles are shown in Figs. 11-13, respectively. We 
note five items of interest in these figures. (1) In 
general, all turbulence models produced very simi- 
lar velocity profiles at each of the stations; the slat 
wake was (as expected) predicted to be too wide 
and too deep, although SARC decreased the wake 
depth slightly. Many other features of the experi- 
mental profiles were captured extremely well by all 
models. ( 2) T he EASMCC overpredicted the mag- 
nitude of u'v* near d/c — 0.01 in Fig. 11; however, 
this overprediction is believed to be due to the use 
of the Laplacian smoother so near to a region of 
high Sij gradient. (3) All models underpredicted 
the peak u'v* near d/c — 0.02 at midflap (Fig. 13). 
(4) Both SARC and EASMCC predicted some lo- 
cal differences in u'v ' levels, but the curvature cor- 
rections overall had relatively minor effects. (5) As 
noted earlier, the EASM did a better job predicting 
the turbulent shear stress just downstream of the 
main element trailing edge [x/c— 0.898); the EASM 
also did a better job predicting the same quantity in 
the region of the upper half of the main wake (near 
d/c — 0.04) at midflap (Fig. 13). It was previously 
thought that poor prediction by the SA model here 
possibly was due to the model’s inability to account 
for curvature effects. However, S ARC’s curvature 
correction did not have much impact. As shown in 
Fig. 14, the nonlinear terms appear to have the most 
influence on the u'v' levels in the midflap area. 

4.4 Global Effects near Maximum Lift 

Clearly, many of the discrepancies in turbulent 
shear stress between previous CFD results 2 and ex- 
periment in the current region of focus were due 
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to use of the Boussinesq eddy viscosity hypothesis. 
Nonetheless, these discrepancies appear to have little 
effect on the velocity profiles in this region: EASM 
yielded generally improved turbulent shear stress 
predictions over the flap, but in spite of this improve- 
ment, the model did no better overall than SA in the 
prediction of mean velocity. Curvature corrections 
in both the one-equation and two-equation models 
had a relatively minor effect for this flow. How- 
ever, the fact that the SARC decreased the depth 
of the slat wake is intriguing and merits more de- 
tailed study. 

One important question remains to be answered. 
Recalling past studies (e.g., Ref. 28) where, com- 
pared with nominally 2-D experiments, CFD over- 
predicted maximum lift, we now ask whether any 
of the turbulence model improvements described in 
this paper improve the global comparisons. Figure 
15 shows results for SA, SARC, and EASM at an- 
gles of attack beyond a = 19 deg. The EASMCC 
was not run for these cases. As in the past stud- 
ies, all three models yielded higher Ct levels and 
a higher stall angle (a = 23 deg) than experiment 
(a = 21 deg). In other words, improved turbulence 
shear stress predictions of EASM had little effect on 
CFD global force results near maximum lift. Curva- 
ture correction terms in SARC also had little effect. 
The reason for CFD’s overprediction of maximum 
lift levels remains unknown. However, results from 
this study support the conclusion from an earlier 
study 1 that turbulence modeling is probably not the 
primary cause for the disagreement. A more likely 
cause is 3-D effects at high angles of attack in the 
wind tunnel; such effects are obviously not modeled 
by 2-D CFD. 

5 Summary 

Several turbulence model improvements, devel- 
oped from unit problem investigations, were incor- 
porated into a study of the flow over the 30P-30N 
multielement airfoil. As a result, turbulent shear 
stress prediction capabilities have been improved, 
and the influence of certain turbulence model com- 
ponents has been established. The K-e form of the 
EASM was shown to be ill suited for use in this flow 
field because it overpredicts turbulence in regions of 
adverse- pressure-gradient wall-bounded flow. This 
problem is inherent in the “standard” form of the 
K-e equations; the K-uj form does not suffer from 
this problem. A new K-uj form of the EASM was 
introduced. Nonlinear terms, present in the EASM, 
were shown to improve the prediction of the turbu- 
lent shear stress behavior behind the trailing edge 
of the main element and near midflap, but these im- 


provements had little effect on the mean flow field. 
Finally, curvature corrections in two different mod- 
els were applied to this flow field. Overall, these cor- 
rections did not have a significant effect in the flap 
region, where the flow field undergoes the greatest 
curvature. Predictions of maximum lift were essen- 
tially unaffected by the turbulence model variations 
studied. 
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Appendix 

This appendix explicitly defines the EASM(AT-u>) 
turbulence model, for readers interested in coding 
it themselves. For details on the derivation of the 
EASM, the reader is referred to earlier papers on 
EASM listed in the references. 

The turbulent stress tensor for EASM is given by 

r ij — — 2i/ { (^Sij 

+ [a 2 a 4 (SikW k j - WikSkj) 

“2a 3 a 4 (SikSkj - ^SkiSik £ij)] ) , (9) 

where 5»j = [{dui/dxj) 4- ( duj/dxi)]/2 and Wij — 
[( dui/dxj ) — {duj/dxi)\f 2. The nonlinear terms are 
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within the brackets []. The component terms are 
used to close the Reynolds- aver aged Navier-Stokes 
equations (see, e.g., Ref. 29). The kinematic eddy 
viscosity is given by 


If d < 0, then 


v? = C*Kt = -Ka i, 


( 10 ) 


with r = 1 /w. Thus, a.\jr is equivalent to — C*. 
The value of ai/r is obtained from the solution to 
the following cubic equation at each point in the flow 
field: 


where 




( 12 ) 


( 277 ^ 27* )2 ( 7l ‘ 2 

+ 27^ 2 ?/ 2 r 2 a2^ 


r = 


7 1 a 1 


(2r7 2 r 2 7*) 2 ' 


root with the lowest real part. 30 Also, the degener- 
ate case when rf* — ► 0 must be avoided. The current 
solution procedure used is as follows: 

If 7j 2 r 2 < 1 x 10" 6 , then 


= 7i Q i n ^ 

\ T J 7 i * 2 -2{W 2 }r 2 al’ 1 J 


Otherwise, define 


2 


_ P 

a = q ~ 3 


6= 27 (2p 3 — 9p<7 + 27r) 

_ ft 2 a 3 
d= J + 27' 


(16) 

(17) 

(18) 


Then, if d > 0 

/ A \ J /3 

<l== V 2 + / ( 19 ) 

/ft N 1/3 

*2 =(-3-^) (20) 

( t ) = ™" (~3 + + |J ' ~3 - 7 - 7 ) (21) 


= cos 1 


2' V 27 


<i = _ o + 2 cos - 



<2_ “3 +2 V I COS 


/27T 0 

v 3 + 3 

/4tt 0 

(t + 3 

(^7) = min(< 1 ,< 2 ,<3) • 


<3 = -f + 2 \H cos 


(22) 

(23) 

(24) 

(25) 

(26) 


In the current implementation, the resulting C* = 
— (ai/r) is limited by C * — max(C7*, 0.0005). The 
nominal level for C * in a zero-pressure-gradient log 
layer is approximately 0.09. 

Other parameters are given by 


17 2 = {S 2 } = SijSji = SijSij 


(13) 

{W 2 } = Wij 


n 2 = 

(14) 

1 


ai = 2 


{W 2 } 


{S 2 } 

\-C7 


02 = 2 (2 - C A ) 

a 3 = 2 (2 - C3) 


Also, 


“4 = [71 ~ 2tS (^f) »7 27 ' 2 ] 


75 = Cl/2 

_ I C 0 , gf 2 ~ Ctl 

71 - 2 Cl + i, <? tl - 1 


(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 


and C c 1 = 1.44, C e2 = 183, C? = 3.4, C\ = 1.8, 
Ci - 0.36, C 3 = 1.25, and C 4 = 0.4. 

The above implementation is exactly the same for 
EASM(A'-w) or EASM(A'-£), except that r = K/e 
for EASM(A'-£). However, the models are different 
with regard to the two-equation model to which they 
are coupled. For EASM(A'-a>), the explicit tensor 
representation for r,j is coupled with the following 
K-u two-equation model: 


DK ^ r „ 9 

- = V-f 0 .Ku+ — 


Duj 

~Dt 


CJ 


= 7—^ - flu + 


d 

dx k 


[( 

[(■ 


OK 


d_ K_ 

dx k 

duj 

dx k 


(36) 

(37) 
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where 


(38) 


p = - nj i)h 2 ‘'‘ v 


and <r k = 2, <t„ = k 2 /[</C^(0 - 7)], « = 0.41, 
7 = 0.575, /? = 0.83, and = 0.0895. Note that 
for 2-D incompressible flows, V = 2z/*?7 2 is exact. In 
the current implementation, V in the /^-equation is 
limited to be less than 20 times the destruction term 
fp'Ku). The function from Wilcox 22 improves 
the performance of the K-u model for 2-D shear 
layers, wakes, and jets and is given by 


u- = 


//5* = 1 when 

1 + 680yg 


1 + 400Xfc 


Xfc < 0 
when Xk > 0 


Xk ~ 


Cl 8K_ du> 
u 3 dxj dxj ' 


(39) 

(40) 

(41) 


where the C 2 term in the formula for \k is necessary 
because w in the current model does not “absorb” 
Cp as in Wilcox’s model. 

The boundary conditions applied at solid walls 
are K w — 0 and uj w = 10(6^) /[/?(An) 2 ], where 
An is the distance to the first cell center away from 
the wall. The boundary condition for ut w is from 
Menter. 5 This boundary condition simulates the an- 
alytical behavior of w near solid walls without the 
need for specifying the solution at interior points. 
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Figure 2. Velocity profiles near trailing edge of main 
element at x/c = 0.85. 




Figure 1. The 30P-30N multielement configuration. 


log(y+) 


Figure 3. Velocity profiles in wall variables at 
x/c = 0.85. 
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Figure 4. Experimentally measured turbulent shear 
stress at two locations near the main element trailing 
edge (dashed curve shifted down by 0.007 units). 



u'v' 


Figure 6. The EASM turbulent shear stress at 
two locations near the main element trailing edge 
(dashed curve shifted down by 0.007 units). 



u'v' 


Figure 5. The SA model turbulent shear stress at 
two locations near the main element trailing edge 
(dashed curve shifted down by 0.007 units). 



u'v 1 


Figure 7. Linear EASM turbulent shear stress at 
two locations near the main element trailing edge 
(dashed curve shifted down by 0.007 units). 
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Figure 8. Velocity profiles near trailing edge of main Figure 10. Velocity profiles near midflap at 
element at x/c = 0.85. x i c = L032 ‘ 



Figure 9. Velocity profiles on flap just downstream Figure 11. Turbulent shear stress profiles near trail- 
of main element at x/c — 0.898. ing ec |g e 0 f ma i n element at x/c — 0.85. 
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uv' 

Figure 12. Turbulent shear stress profiles on flap 
just downstream of main element at x/c = 0.898. 



u'v* 


Figure 14. Effect of nonlinear terms on turbulent 
shear stress profiles at x/c = 1.032. 



uV 

Figure 13. Turbulent shear stress profiles near mid- 
flap at x/c = 1.032. 



a , deg 

Figure 15. Lift coefficients near maximum lift. 
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